The role of electronic excitation in cold atom-ion chemistry 
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The role of electronic excitation in charge exchange chemical reactions between ultracold Ca 
atoms and Ba + ions, confined in a hybrid trap, is studied. This prototypical system is energetically 
precluded from reacting in its ground state, allowing a particularly simple interpretation of the 
influence of electronic excitation. It is found that while electronic excitation of the ion can critically 
influence the chemical reaction rate, electronic excitation of the neutral atom is less important. It 
is also experimentally demonstrated that with the correct choice of the atom-ion pair, it is possible 
to mitigate the unwanted effects of these chemical reactions in ultracold atom-ion environments, 
marking an important step towards the next generation of hybrid devices. 

PACS numbers: 




Since the inception of laser-cooling, the primary focus 
of atomic physics has been the development of techniques 
for the production and study of ultracold matter - an 
endeavor that, at its core, is centered on gaining full con- 
trol over matter at the quantum level. This work has 
been extremely successful, enabling many long-sought- 
after goals, such as quantum degenerate gases^P, 
quantum simulation [3(, quantum information 
precision measurement of fundamental physics 
gain this control, however, most ultracold matter produc- 
tion techniques rely on the scattering of a large number of 
photons from the system under study, potentially leading 
to a large degree of electronic excitation. Thus, the atom 
or molecule being cooled can be in a somewhat peculiar 
state: its external motion may be characterized by a tem- 
perature close to absolute zero, but its internal electronic 
degree of freedom may be described by a temperature ap- 
proaching infinity. While this non-thermal distribution 
of electronic states has some notable important conse- 
quences for ultracold atoms Ml, most of its effects, such 
as photochemical reactions (si— 1 1 01] - are largely ignored as 
they occur at rates that have relatively little effect on the 
system. However, as the field now moves towards pro- 
ducing more complex systems at ultracold temperatures, 
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e.g. molecules 1 1 II ] and hybrid systems 
feet of these light-assisted processes must be reevaluated 
as their rates may be much larger due to, among other 
things, an increased density of accessible product states 
and longer range interactions. Thus, there is presently a 
need to better understand the role of electronic excitation 
in chemical reactions at ultracold temperatures to enable 
the next generation of atomic physics experiments. 

Interestingly, the rapidly emerging field of hybrid 
atom-ion systems offers a unique opportunity to study 
ultracold chemical reactions |14l4l7| . Like the all-neutral 
systems of traditional atomic physics, the two trapped 
species can collide at short-enough range for chemical 
reactions to proceed; but, unlike all- neutral systems it 
is possible to maintain a product of the chemical reac- 



tion in the trap, since ion trap depths are large rela- 
tive to the kinetic energy gained in most chemical re- 
actions. Already, such hybrid systems have been used 
to measure several important ion-neutral chemical reac- 
tions culminating in the recent observation of 



molecular ion reaction products [17j . Here, we utilize a 
hybrid atom-ion MOTION trap system [l6, 17 1 to study 
reactions between ultracold 40 Ca atoms and 138 Ba + ions. 
Unlike all previous studies, the only energetically acces- 
sible chemical reaction pathways for this system require 
the atom, the ion, or both to be electronically excited. 
As this reaction takes place between a closed-shell atom 
and open-shell ion, the results of this system are simpler 
to interpret and serve as a general prototype for draw- 
ing conclusions about more complicated systems. Fur- 
ther, by utilizing the richer electronic structure of 40 Ca, 
as compared to neutral atoms used in other hybrid sys- 
tems, we are able to vary the population in the neutral 
atom electronic states over a wider range than in previous 
experiments and present the first detailed experimental 
investigation of the importance of neutral atom excita- 
tion in ion-neutral chemistry. 

In what follows, we briefly describe our experimental 
system, discuss how optical coherence and trap dynam- 
ics can critically affect the observed atom-ion chemical 
reaction rate when electronic excitation is involved, and 
present experimental results detailing the observed chem- 
ical reaction pathways and the average radiative associ- 
ation branching ratio for formation of CaBa + molecu- 
lar ions. We also present relativistic ah initio molecu- 
lar structure calculations, which further illuminate the 
chemical reaction pathways as well as provide guidance 
for future experiments. We conclude with a demonstra- 
tion that it is possible to avoid the large chemical reaction 
rates seen in other hybrid atom-ion systems, with rea- 
sonable experimental parameters, marking a crucial step 
towards the next generation of hybrid atom-ion devices. 

The apparatus used in this work is similar to that de- 
scribed in Ref. [l6| . Briefly, 138 Ba + ions are loaded into 
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FIG. 1: a) Ba + and Ca relevant energies levels and laser cooling schemes, b) Fluorescent images of Ca MOT and Ba + ion cloud, 
c) and d) Molecular potentials of the f2=|, | excited states of the BaCa + molecule, respectively, as a function of interatomic 
separation R. The bottom panels show higher detail for potentials near 20 x 10 3 cm -1 above the Ca + ( 2 S)+Ba( 1 D) dissociation 
limit. Entrance and exit channels, as defined in the text, dissociate to green and red potentials, respectively. 



a linear radio-frequenc y q uadrupole trap via laser abla- 
tion of a BaCl2 target [18[ and cooled with colinear laser 
beams. In the same region, ultracold Ca atoms are pro- 
duced and trapped in a magneto-optical trap (MOT). 
(See Fig. [lja) for the respective ion and atom laser cool- 
ing schemes.) Two nearly-orthogonal cameras allow for 
a 3-dimensional reconstruction of both the Ca MOT and 
the laser-cooled ion cloud, and are used to quantify the 
degree of overlap between the Ca atoms and 138 Ba + ions, 
as shown in Fig. Eh). Using the same approach as 
Ref . [l6| , the laser-induced 493 nm fluorescence monitors 
the 138 Ba + ion population, which is observed to decay 
in the presence of the Ca MOT due to charge exchange 
chemical reactions. Special care must be taken to account 
for the time evolution of the degree of overlap, which re- 
quires a non-analytical numerical fitting technique, de- 
scribed in Appendix A, to account for the deviation from 
simple exponential decay. 

In this experiment, the atom- ion collision energies are 
primarily set by the ion micromotion since the temper- 
ature of the Ca MOT (T ~ 4 mK) and the Doppler- 
cooling temperature limit of the Ba + ions (T ~ 0.5 mK) 
are both lower than the typical, position dependent, en- 
ergy of micromotion in the ion cloud ~10 mK-10 K). 
At these low collision energies, chemical reaction dynam- 
ics are primarily determined by the relative energies of 
the entrance and exit reaction channels, shown in black 
and red, respectively, in Fig. [Hc,d). Because the ioniza- 
tion potential of Ca (6.1 eV) is significantly larger than 
that of Ba (5.2 eV), charge exchange reactions between 
the Ca 4 1 S- and Ba + 6 2 S- or 5 2 D- electronic states are 
energetically precluded. As a result, chemical reaction 
between Ba + and Ca can only occur if the atom, the ion, 
or both are electronically excited. Of the ten energeti- 
cally allowed entrance channels populated in this experi- 



ment, the six involving the short-lived Ca P-states do not 
contribute to the observed reaction rate. Because the Ca 
4 1 S-, 4*P- and 5 1 P-states have very different atomic po- 
larizabilties [l9| (163.0 a.u., 55.3 a.u., and >1000 a.u., re- 
spectively), during the collision the strong monopolc field 
of the ion alters the atomic transition frequency such that 
the Ca MOT lasers are shifted out of resonance for atom- 
ion separations of many hundred Bohr radii, a . Thus, 
though a collision leading to chemical reaction may begin 
with a Ca atom in a P-state, by the time the atom and ion 
are close enough together to react (< 10 a ), the Ca atom 
has radiatively decayed from the P-state. We note similar 
effects are discussed in Ref. [3, 21 1 , were observed in 
Ref. [lfj ] , and are confirmed by the data presented here. 
This effect appears generic and it is likely that short-lived 
excitations of ultracold atoms can be ignored in future 
experiments involving reactions of ultracold atoms and 
ions. (This effect is particularly beneficial to techniques 
for sympathetic cooling of molecular ions by ultracold 
atoms [221 ] . as it significantly relaxes the experimental 
constraints.) Interestingly, a similar argument cannot be 
made for short-lived, excited states of ions. Because ex- 
citation of the ion only changes the long-range dispersion 
coefficient, e.g. Cg, the ion cooling lasers are not shifted 
out of resonance until relatively small atom-ion separa- 
tions. Thus, while these chemical reactions are some- 
what suppressed, an effect that must be accounted for 
when evaluating rate constants involving excited states 
susceptible to spontaneous decay, it is nonetheless possi- 
ble for them to occur between atoms and excited state 
ions. 

With these considerations in hand, it is clear that the 
important reaction pathways proceed from the energet- 
ically allowed entrance channels involving 4 1 S- or 3^- 
state Ca atoms and 6 2 S-, 6 2 P-, or 5 2 D-state Ba + ions, 
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FIG. 2: Panel (a): experimentally determined excitation frac- 
tion in Ba + versus theoretical predictions of a density matrix 
(•) and rate equation (o) treatment. Panel (b): experimen- 
tally determined excitation fraction in Ca versus predictions 
of a 4-level rate equation model. In both panels the solid line 
has a unit slope and represents perfect agreement between 
experiment and theory. 



highlighted in green in Fig.QJc-d). While the total charge 
exchange rate constant for each of the four distinct path- 
ways depends sensitively on the details of the CaBa + 
molecular structure, it can nevertheless be measured in a 
straightforward way by simply varying the excited state 
populations of the atoms and ions while recording the to- 
tal system reaction rate. Experimentally, this variation is 
accomplished by changing the intensities of the atom and 
ion cooling and repumping lasers. From the known ex- 
perimental parameters it is then possible to calculate the 
relative population of each atom and ion electronic state, 
and extract the contribution of each reaction channel to 
the total measured reaction rate. 

The calculation of the Ba + electronic state populations 
is significantly complicated by optical coherence effects 
and the quadrupole magnetic field of the MOTION trap, 
whose direction and strength vary in a non-trivial way. 
Therefore, we determine the Ba + electronic state pop- 
ulations by solving the steady-state quantum Liouville 
equation for an eight-level optical Bloch Hamiltonian at 
each position in the ion cloud and spatially averaging 
the result, as described in Appendix B. Similarly, a first 
principles calculation of the electronic state populations 
of atoms in a MOT environment is a well-known open 
problem and various phcnomenological models have been 
proposed |23j,[24|. Thus, to verify the calculation of Ba+ 
populations and determine the appropriate model for Ca 
populations, we experimentally measure the fraction of 



population in both the Ba + 6 2 Pi/ 2 state, and the Ca 
4 1 P-state, Ppp, by comparing the observed fluorescence 
rates of each species to the total number of scatters as 
measured by a channeltron ion detector and absorption 
imaging, respectively. The results of the Ba + popula- 
tion measurement are summarized in Fig. [Ha), where 
the predictions of a simple rate equation model are in- 
cluded for comparison. The poor agreement of the p^f 
predicted by the rate equation model highlights the need 
for the density matrix treatment, as its use can lead to 
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FIG. 3: Measured rate constant as a function of Ba + popu- 
lation fractions. Experimental result in black. Best fit to the 
four entrance channel model shown as open red circles (o). 




FIG. 4: Measured rate constant as a function of Ca population 
fractions. Experimental result in black. Best fit to the four 
entrance channel model shown as open red circles(o). 



significant underestimate of the charge exchange reaction 
rate constant in three-level systems. The results of the 
Ca population measurement arc shown in Fig. EJb) along 
with the prediction of a simple four-level rate equation 
model. The good agreement with the model is likely 
a consequence of the lack of hypcrfine structure in the 
40 Ca system; however, in the absence of a first principles 
model, it could simply be fortuitous. Using this result, 
we estimate the fraction of Ca in the 3 1 D-state, p^, by 
finding equilibrium between the ^P— >3 1 D spontaneous 
decay rate and the combination of the ballistic exit rate 
from the MOT and laser repumping rate (lo| . 

With these tools in place, we measure the reaction 
rate constant as a function of the Ba + and Ca excited 
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state fraction, as shown in Figs. [3] and [4] respectively. 
Care must be taken in interpreting the measurements 
presented in these graphs due to the correlations among 
the fractional populations as the laser intensities are var- 
ied. As a function of Ba + excitation, the most obvi- 
ous feature is a linear trend of increasing reactivity with 
increasing Pp p a . In addition, the reaction rate constant 
diminishes at both extremes of the Ba + S-state, pf s a , 
and D-state, p^, population fraction. Taken together, 
these features indicate that the primary reaction chan- 
nel involves the Ba + ion P-statc. Also shown in Fig. G2 
as open triangles, are several measurements at low exci- 
tation fractions (< 1%). While these points appear to 
qualitatively agree with the rest of the data, we report 
data analysis both with and without inclusion of these 
points in order to avoid possible systematic inaccuracies 
due to higher secular energies. 

The measurement of the rate constant over the experi- 
mentally accessible Ca excitation range, Fig. HI shows no 
detectable variation for fixed Ba + trapping parameters. 
With the aforementioned assumption that Ca atoms in 
the 4 1 P-state radiatively decay to the 4 1 S-state by the 
time they are close enough to an ion to chemically react 
and with negligible 5 1 P-state population for our operat- 
ing parameters, this data indicates that the Ca S-state 
and D-state exhibit similar reaction rates with Ba + ions. 
However, since the four-level rate equation model has 
only been verified to accurately predict the fraction of 
Ca atoms in the P-state (experimental verification of the 
D-state populations in not presently possible in our ap- 
paratus) , the possibility remains for systematic errors as- 
sociated with an inaccurate calculation of the Ca D-state 
population. 

To extract channel-specific reaction rate constants 
from these data, we quantify the total reaction rate con- 
stant in terms of the excited state fractions of both Ba + 
and Ca as: k = p^ p a p^ a k ps + pf s a p^k sd + pftpftkM + 
PppPdd^pd- By accounting for the conservation of prob- 
ability (p ss = 1 — p pp — pdd) and the fact that Ca atoms 
in the P-state radiatively decay to the S-state before 
the atom and ion are close enough to react {p* s f a = 
pg a + p%°), the totality of the data in Fig. Hand Fig. H 
was subject to a multi-parameter fit along the manifold 
{Pp P a ,Pdd>Pdd} t° determine the contribution of each 
reaction pathway to the total rate constant. The re- 
action rate constant of the primary entrance channel, 
Ba+(6 2 Pi /2 ) + Ca(4 1 S), was measured to be kp S — 
4.2(1.9) x 10~ 10 cm 3 Hz. The constraints for the other 
entrance channels vary depending on whether the low ex- 
citation fraction Ba + (6 2 P!/ 2 ) data is included in the fits. 
If the small excitation population points are included in 
the analysis, the rate constant for the only channel out of 
the Ba + ground state, Ba + (6 2 Si/ 2 ) + Ca(3 1 D), is poorly 
constrained with a best fit of k s d ~ 2.0 x 10~ 10 cm 3 
Hz. If the points are not included, the best fit for k s d 



becomes consistent with zero within a large standard 
error of « 3.0 x 10~ 10 cm 3 Hz. Because the remain- 
ing two entrance channels, Ba + (6 2 Pi/2) + Ca(3 1 D) and 
Ba + (5 2 D 3 / 2 ) + Ca(2 1 D), require excitation of both the 
atom and ion our experiment is relatively insensitive 
to them and we can only set upper limits on them of 
k p d < 2 x 10~ 9 cm 3 Hz and k dd < 7 x 10~ 10 cm 3 Hz. 
This fitting method yields a surface in dimensions of the 
various excitation fractions which was traversed in a non- 
trivial manner under variation of laser intensities, and 
thus it is not possible to plot a line of best fit. There- 
fore to show the quality of agreement between the best 
fit model and the results, we plot a data point gener- 
ated from the model corresponding to each experimental 
value (open circle, red) , with error bars including the full 
uncertainty associated with the upper bounds of the dou- 
bly excited channel. The relatively uncertain agreement 
for the measurements taken at higher Ca D-state popu- 
lation is indicative of either a failure in our estimation 
of the D-state fraction or that the upper bounds are too 
conservative. 

To illuminate the charge exchange mechanisms respon- 
sible for the individual reaction pathways, we performed 
ab initio calculations of the 0=1/2, 3/2 potentials us- 
ing a relativistic multi-reference restricted active space 
configuration- interaction method [25|. Spin-orbit effects 
are large in CaBa + and a relativistic calculation is re- 
quired as described in Appendix C. The resulting Born- 
Oppenheimer potentials for the ground and first-excited 
state of CaBa + are shown in Fig. Qic-d), an d assigned 
by their atomic dissociation limit. The bottom panels of 
Fig- Q3c-d) show a higher resolution view of the closely 
spaced potentials above 20 x 10 3 cm -1 . Both panels re- 
veal that there are strong interactions and avoided cross- 
ings between neighboring excited potentials. Potentials 
dissociating to energies > 27 x 10 3 cm -1 are not shown, 
as the exceedingly high level density makes it impossible 
to resolve the avoided crossings with our calculation. 

In principle, a coupled-channels calculation that in- 
cludes the mixings between these potentials due to ra- 
dial non-adiabatic coupling and spontaneous emission, 
like that performed in Ref. [3], can be used to calcu- 
late the non-radiative, radiative, and radiative associ- 
ation charge exchange reaction rate constants for each 
entrance channel; however, due to the high level den- 
sity in the excited CaBa + molecule, this calculation is 
extremely technically demanding and falls outside the 
scope of this work. Nonetheless, the basic mechanisms 
of each charge exchange pathway can be inferred from 
the structure of its entrance channel molecular potential. 
For example, the experimentally determined primary en- 
trance channel, Ca(4 1 S) + Ba + (6 2 P!/ 2 ), exhibits strong 
interactions with several potentials as well as a narrowly- 
avoided crossing with a molecular potential dissociating 
to the Ca+(4 2 S) + Ba(6 3 P) exit channel. The second 
experimentally addressed entrance channel, Ca(3 1 D) + 
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Ba + (6 2 Si/2), has several strong avoided crossings with 
the Ca + (3 2 D) + Ba(6 1 S) exit channel. These qualities 
imply that non-radiative charge exchange is most likely 
responsible for the observed reaction rate. By a similar 
analysis, it is expected that the Ca(3 1 D) + Ba + (5 2 D) en- 
trance channel would be relatively inefficient for charge 
exchange reactions since low energy scattering events 
do not experience an avoid crossing nor do they cou- 
ple well to states accessible via spontaneous emission. 
The remaining Ca(3 1 D) + Ba + (6 2 Pi/2) entrance chan- 
nel resides in a host of densely packed levels, which both 
prevents our ab initio calculation from providing reli- 
able results and suggests that non-radiative and radia- 
tive chemical reactions could proceed very quickly from 
these states, a result that is consistent with the limit set 
for this channel by our measurement. 

Despite the fact that the rate constants measured here 
are large by atomic physics standards, roughly a third of 
the Langevin rate, with proper engineering of the atomic 
and ionic excited state populations their effects can be 
mitigated. By using large repumping laser powers to 
maintain a negligible Ca D-state fraction and using the 
minimum ion laser cooling power necessary to achieve 
Dopplcr-limited cooling, we have observed Ba + ion life- 
times in excess of several minutes in the presence of the 
MOT. Finally, using the method described in Ref. flij ]. 
we measure the molecular ion product branching ratio 
to be ~ 0.14. Future measurements, using time-of-flight 
detection as in Ref. [26j , are planned to further constrain 
the branching ratio. 

In summary, we have experimentally and theoretically 
studied chemical reactions in a hybrid atom-ion system, 
where only electronically excited states are permitted to 
react. The observed reaction dynamics are qualitatively 
different than what is observed in systems where ground- 
state reactions dominate fl6| . The measured reaction 
rate constant is found to be an appreciable fraction of 
the Langevin rate and most-likely due to non-radiative 
charge transfer occurring at narrowly-avoided crossings 
in at least two of the entrance channels. We have also 
demonstrated that while special care must be taken to 
account for e.g. optical coherence and trap dynamics 
when understanding the degree of electronic excitation, 
it is possible to engineer experiments where the cold reac- 
tive species can co-exist for several minutes. From these 
results it is clear that the combination of ultracold atoms 
with high ionization potential and ultracold ions with low 
electron affinities, along with careful control of electronic 
excitation should enable a new generation of hybrid-atom 
ion devices, capable of the long coherence times needed 
for proposed hybrid atom-ion systems 
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AFOSR grant No. FA 9550-11-1-0243. 
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Appendix A: Time-varying overlap factor 

As the Ca MOT exhibits non-uniform particle densi- 
ties on length scales comparable to the size of laser-cooled 
ion cloud, it is necessary to carefully consider the effec- 
tive density of rcactants when evaluating the reaction 
rate constant. Further, because the chemical reaction 
leads to loss of ions, the ion cloud size shrinks during the 
reaction process, resulting in a time- variation of this ef- 
fective density. This effect can cause various systematic 
shifts in the measurement of the rate constant if not taken 
into account. For example, if the ion and atom clouds are 
nonconcentric, as the ion cloud shrinks the observed reac- 
tion rate slows as a result of the reduced reactant density. 
Thus, to properly measure the reaction rate constant, we 
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explicitly include the time evolution of the effective reac- 
tant density via a term that quantifies the spatial over- 
lap of the reactants, defined as </>(£) = J /5 Ca (r)p/(r)dr, 
where p denotes a peak normalized density and p denotes 
an integral normalized density. With this definition, the 
evolution of the number of Ba + , N, ions is given as: 



dN 
~dt 



-kp<j>(t)N 



(1) 



where we take the product of the reaction rate constant, 
k, to the peak neutral density, p, to be the unbiased 
reaction rate constant Tcex- 

In principle, <p(t) can be found numerically from mea- 
surements of the cloud widths and relative positions via 
dual camera fluorescence imaging at various times. How- 
ever, due to disparate levels of brightness, the MOT and 
the ions cannot be simultaneously imaged with our cur- 
rent image filtering systems. Therefore, in practice, we 
measure the relative position of the two clouds and their 
initial density profiles. With the assumption that the Ca 
MOT has a roughly spherical gaussian density distribu- 
tion, while the ion cloud has an ellipsoidal (cigar-shaped) 
gaussian density distribution, we can then uniquely de- 
termine numerically the time evolution of the overlap 
factor (in units of Tq^), as in Fig. [5] Thus, for the de- 
termination of Tcex for each decay measurement, a nu- 
merical function describing the overlap factor is uniquely 
determined from the images so that the only physical 
fitting parameter is Tcex- That is, each experimental 
decay curve is effectively fit to ^ = — FcEX</>(rcExi)-^ 
with an uniquely determined </>. In order to use numer- 
ical fitting algorithms, we take the solution to the dif- 
ferential equation to be N(t) ~ exp(— rcEX'MrcExt))- 
This amounts to an assumption that the overlap factor 
is slowly varying on the timescale of the reaction and for 
the data reported here leads to systematic error < 30%, 
which is well within the reported experimental error. 



Appendix B: Calculation of Ba + fractional electronic 
state populations 

Because the trapped Ba + ions are laser cooled via a 
closed three- level electronic system, see Fig. [TJa), it is 
necessary to account for optical coherence effects to ac- 
curately calculate the excited state population fractions 
realized in the experiment. This is accomplished by solv- 
ing for the steady-state solution to the quantum Liou- 
ville equation for an eight-level optical Bloch Hamilto- 
nian. This calculation is significantly complicated for the 
MOTION trap due to the presence of a large quadrupole 
magnetic field necessary for the operation of the Ca MOT 
as its direction and strength vary in a non-trivial way. We 
briefly outline this calculation in this appendix. 

Working under the electric dipole approximation in the 
interaction picture, the eight bare \{LS},J,Mj) states 




FIG. 5: The dashed line shows the numerical calculation of 
the overlap factor as a function of charge exchange lifetimes 
expired in the case of perfect alignment of relative MOT to 
ion cloud position, with typical experimental waists and as- 
pect ratio. The other three lines are randomly selected nu- 
merical calculations from actual experimental data to show 
how imperfect alignment changes the curve. 



have one of three energies, E$ = 0,Ep = S,Ed = A, 
where S and A are the one- and two-photon detunings, 
respectively. To simplify the geometry of the problem, 
we define the quantization axis to be parallel to the di- 
rection of the linearly polarized laser used to drive the 
Ba + 5 2 Pi/2 <— 5 2 S cooling transition, which is approxi- 
mately vertical in the laboratory frame. The laser used 
for population repumping on the 5 2 Pi/2 4— 5 2 D 3 / 2 tran- 
sition is collinear with this cooling laser and is linearly 
polarized at an angle 9 with respect to the z-axis. As a 
result, the Rabi frequencies for transitions between the 
different \{LS}J,Mj) states of the Ba+ ion are given as: 



Sl L ,, L (J',M'j,J,Mj,E ltq 
J' 1 J 



(2) 



( - 1)J_Mj { M'j q MjjiL'iJ'WdimJ^Jh 

where the index L = 2 S, 2 Pi/2, 2 D 3 / 2 and £\ t o = £ z 
and £i,±i = ±-^(£ x T l £y) are the spherical irreducible 
tensor representation of the electric fields of the laser res- 
onant with the L'j, <— Lj transition - with the defined 

geometry, £ x = £ y = 0, £ z = J Mes. for the cooling laser 



and is £ x 



Me? sm6, £ v = 0,£ z 



Mell C o S for 



the repumping laser, where Ipl is the intensity of the 
laser driving the transition. The angle 9 is experimen- 
tally chosen to mitigate the effects of coherent population 
trapping and, in our experiment, is typically 9 as 7r/3. 

Likewise, the effect of the spatially varying magnetic 
field is calculated by including in the Hamiltonian the 
diagonal and off-diagonal Zeeman matrix elements: 



(J',M'j\-p-B\J,Mj) = 

eh , „ / J+Bi _i J-Bi i 



(3) 



2m, 



g{J) 
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where B lfi = B z . B lt ±i = ^-±=(B X ± iBy), and g(J) = 
(3J(J + 1) + S{S+ 1) - L(i + 1))/(2J(J+1)). Over the 
spatial extent of the Ba + ion cloud the magnetic field 
is well approximated as B = fr(f x + |y + zz), with 
4^ = 0.7 T/m for our experimental parameters. Here we 
ignore off-diagonal Zeeman couplings between different 
electronic states, e.g. matrix elements between P3/2 and 
Pi/2 states. 

Finally, the effects of spontaneous emission and the re- 
sulting dccohcrence are included via the Lindblad super- 
operator and result in a relaxation matrix, T r , whose di- 
agonal elements account for the decay and growth of elec- 
tronic state populations at the respective spontaneous 
emission rate and whose off-diagonal elements account 
for the concomitant decay in coherence at half of the 
spontaneous emission rate. 

With all of these terms taken together, the optical 
Bloch Hamiltonian, 

ih% = [H(x, Ips,I PD ,9,S,A),p}+ ihT r , (4) 
at 

is then solved numerically for given experimental condi- 
tions, i.e. a given set of {Ips, Ipd, S, A}, to find the 
steady state Ba + populations at each location in the ion 
cloud. A spatial average is then performed over these 
values to give the average electronic state populations of 
the Ba + ions for the experimental conditions. For the 
conditions studied in this work, this technique gives ex- 
cellent agreement with measured values of Ba + electronic 
populations, as evidence in Fig. HJa). 

Appendix C: ab initio calculation 

This RMR-RAS-CI method is based on partitioning 
the occupied and unoccupied molecular orbitals into sub- 



sets, which define the construction of the molecular wave 
function. These symmetry subsets are: core orbitals, 
which do not participate in the CI procedure; valence 
orbitals, which are occupied and can have single, double, 
and triple excitations; and virtual or unoccupied orbitals, 
which have double excitations. Though the inactive core 
orbitals do not contribute to correlation effects, they cre- 
ate, a so-called, "self-consistent-field sea" in which the 
other electrons move. The valence and virtual orbitals 
are correlated in the CI expansion. Since the dimension 
of the CI molecular wave function increases dramatically 
with the size of the active orbital space, it is necessary 
to apply restrictive measures. 



A non-orthogonal basis set is constructed from nu- 
merical Dirac-Fock atomic orbitals as well as relativis- 
tic Sturmian functions. A symmetric re-cxpansion of 
atomic orbitals from one atomic center to another simpli- 
fies the calculation of many-center integrals. At large in- 
teratomic separations the molecular wave function has a 
pure atomic form that appropriately describes the molec- 
ular dissociation limit. The ls 2 2s 2 2p 6 ..Ad 10 closed shells 
of barium and ls 2 2s 2 2p 6 closed shells of calcium atom 
arc included in the core. An R-dcpcndcnt all-clcctron 
core potential is calculated exactly and included in the 
Hamiltonian. The core-valence basis set is constructed 
from Dirac-Fock core 5s and 5p orbitals for Ba and 3s 
and 2>p orbitals for Ca. Valence Dirac-Fock orbitals are 
6s, 5d, and 6p for Ba and 4s, 3c?, and 4p for Ca, while 
Sturm virtual 7s, 7p, 6d, 8s 8p and 7c? orbitals for Ba 
and 5s, 5p, 4c?, 6s, 6p, and 5c? orbitals for Ca are also 
included. 



